MASTR-seq: Multiplexed Analysis of Short Tandem Repeats with sequencing.

More than 60 human disorders have been linked to unstable expansion of short tandem repeat (STR) tracts. STR length and the extent of DNA methylation is linked to disease pathology and can be mosaic in a cell type-specific manner in several repeat expansion disorders. Mosaic phenomenon have been difficult to study to date due to technical bias intrinsic to repeat sequences and the need for multi-modal measurements at single-allele resolution. Nanopore long-read sequencing accurately measures STR length and DNA methylation in the same single molecule but is cost prohibitive for studies assessing a target locus across multiple experimental conditions or patient samples. Here, we describe MASTR-seq, M ultiplexed A nalysis of S hort T andem R epeats, for cost-effective, high-throughput, accurate, multi-modal measurements of DNA methylation and STR genotype at single-allele resolution. MASTR-seq couples long-read sequencing, Cas9-mediated target enrichment, and PCR-free multiplexed barcoding to achieve a >ten-fold increase in on-target read mapping for 8-12 pooled samples in a single MinION flow cell. We provide a detailed experimental protocol and computational tools and present evidence that MASTR-seq quantifies tract length and DNA methylation status for CGG and CAG STR loci in normal-length and mutation-length human cell lines. The MASTR-seq protocol takes approximately eight days for experiments and one additional day for data processing and analyses. Key points: We provide a protocol for MASTR-seq: M ultiplexed A nalysis of S hort T andem R epeats using Cas9-mediated target enrichment and PCR-free, multiplexed nanopore sequencing. MASTR-seq achieves a >10-fold increase in on-target read proportion for highly repetitive, technically inaccessible regions of the genome relevant for human health and disease.MASTR-seq allows for high-throughput, efficient, accurate, and cost-effective measurement of STR length and DNA methylation in the same single allele for up to 8-12 samples in parallel in one Nanopore MinION flow cell.

Post-transcriptional splicing can occur in a slow-moving zone around the gene.

Splicing is the stepwise molecular process by which introns are removed from pre-mRNA and exons are joined together to form mature mRNA sequences. The ordering and spatial distribution of these steps remain controversial, with opposing models suggesting splicing occurs either during or after transcription. We used single-molecule RNA FISH, expansion microscopy, and live-cell imaging to reveal the spatiotemporal distribution of nascent transcripts in mammalian cells. At super-resolution levels, we found that pre-mRNA formed clouds around the transcription site. These clouds indicate the existence of a transcription-site-proximal zone through which RNA move more slowly than in the nucleoplasm. Full-length pre-mRNA undergo continuous splicing as they move through this zone following transcription, suggesting a model in which splicing can occur post-transcriptionally but still within the proximity of the transcription site, thus seeming co-transcriptional by most assays. These results may unify conflicting reports of co-transcriptional versus post-transcriptional splicing.

A multi-looping chromatin signature predicts dysregulated gene expression in neurons with familial Alzheimer's disease mutations.

Mammalian genomes fold into tens of thousands of long-range loops, but their functional role and physiologic relevance remain poorly understood. Here, using human post-mitotic neurons with rare familial Alzheimer's disease (FAD) mutations, we identify hundreds of reproducibly dysregulated genes and thousands of miswired loops prior to amyloid accumulation and tau phosphorylation. Single loops do not predict expression changes; however, the severity and direction of change in mRNA levels and single-cell burst frequency strongly correlate with the number of FAD-gained or -lost promoter-enhancer loops. Classic architectural proteins CTCF and cohesin do not change occupancy in FAD-mutant neurons. Instead, we unexpectedly find TAATTA motifs amenable to binding by DLX homeodomain transcription factors and changing noncoding RNAPolII signal at FAD-dynamic promoter-enhancer loops. DLX1/5/6 mRNA levels are strongly upregulated in FAD-mutant neurons coincident with a shift in excitatory-to-inhibitory gene expression and miswiring of multi-loops connecting enhancers to neural subtype genes. DLX1 overexpression is sufficient for loop miswiring in wildtype neurons, including lost and gained loops at enhancers with tandem TAATTA arrays and singular TAATTA motifs, respectively. Our data uncover a genome structure-function relationship between multi-loop miswiring and dysregulated excitatory and inhibitory transcriptional programs during lineage commitment of human neurons homozygously-engineered with rare FAD mutations.

Cell type-specific loops linked to RNA polymerase II elongation in human neural differentiation.

DNA is folded into higher-order structures that shape and are shaped by genome function. The role for long-range loops in the establishment of new gene expression patterns during cell fate transitions remains poorly understood. Here, we investigate the link between cell-specific loops and RNA polymerase II (RNAPolII) during neural lineage commitment. We find thousands of loops decommissioned or gained de novo upon differentiation of human induced pluripotent stem cells (hiPSCs) to neural progenitors (NPCs) and post-mitotic neurons. During hiPSC-to-NPC and NPC-to-neuron transitions, genes changing from RNAPolII initiation to elongation are >4-fold more likely to anchor cell-specific loops than repressed genes. Elongated genes exhibit significant mRNA upregulation when connected in cell-specific promoter-enhancer loops but not invariant promoter-enhancer loops, promoter-promoter loops, or unlooped. Genes transitioning from repression to RNAPolII initiation exhibit slight mRNA increase independent of loop status. Our data link cell-specific loops and robust RNAPolII-mediated elongation during neural cell fate transitions.

Spatially coordinated heterochromatinization of long synaptic genes in fragile X syndrome.

Short tandem repeat (STR) instability causes transcriptional silencing in several repeat expansion disorders. In fragile X syndrome (FXS), mutation-length expansion of a CGG STR represses FMR1 via local DNA methylation. Here, we find megabase-scale H3K9me3 domains on autosomes and encompassing FMR1 on the X chromosome in FXS patient-derived iPSCs, iPSC-derived neural progenitors, EBV-transformed lymphoblasts, and brain tissue with mutation-length CGG expansion. H3K9me3 domains connect via inter-chromosomal interactions and demarcate severe misfolding of TADs and loops. They harbor long synaptic genes replicating at the end of S phase, replication-stress-induced double-strand breaks, and STRs prone to stepwise somatic instability. CRISPR engineering of the mutation-length CGG to premutation length reverses H3K9me3 on the X chromosome and multiple autosomes, refolds TADs, and restores gene expression. H3K9me3 domains can also arise in normal-length iPSCs created with perturbations linked to genome instability, suggesting their relevance beyond FXS. Our results reveal Mb-scale heterochromatinization and trans interactions among loci susceptible to instability.

Polygenic burden of short tandem repeat expansions promote risk for Alzheimer's disease.

Studies of the genetics of Alzheimer's disease (AD) have largely focused on single nucleotide variants and short insertions/deletions. However, most of the disease heritability has yet to be uncovered, suggesting that there is substantial genetic risk conferred by other forms of genetic variation. There are over one million short tandem repeats (STRs) in the genome, and their link to AD risk has not been assessed. As pathogenic expansions of STR cause over 30 neurologic diseases, it is important to ascertain whether STRs may also be implicated in AD risk. Here, we genotyped 321,742 polymorphic STR tracts genome-wide using PCR-free whole genome sequencing data from 2,981 individuals (1,489 AD case and 1,492 control individuals). We implemented an approach to identify STR expansions as STRs with tract lengths that are outliers from the population. We then tested for differences in aggregate burden of expansions in case versus control individuals. AD patients had a 1.19-fold increase of STR expansions compared to healthy elderly controls (p=8.27×10-3, two-sided Mann Whitney test). Individuals carrying > 30 STR expansions had 3.62-fold higher odds of having AD and had more severe AD neuropathology. AD STR expansions were highly enriched within active promoters in post-mortem hippocampal brain tissues and particularly within SINE-VNTR-Alu (SVA) retrotransposons. Together, these results demonstrate that expanded STRs within active promoter regions of the genome promote risk of AD.

Schizophrenia-associated somatic copy-number variants from 12,834 cases reveal recurrent NRXN1 and ABCB11 disruptions.

While germline copy-number variants (CNVs) contribute to schizophrenia (SCZ) risk, the contribution of somatic CNVs (sCNVs)-present in some but not all cells-remains unknown. We identified sCNVs using blood-derived genotype arrays from 12,834 SCZ cases and 11,648 controls, filtering sCNVs at loci recurrently mutated in clonal blood disorders. Likely early-developmental sCNVs were more common in cases (0.91%) than controls (0.51%, p = 2.68e-4), with recurrent somatic deletions of exons 1-5 of the NRXN1 gene in five SCZ cases. Hi-C maps revealed ectopic, allele-specific loops forming between a potential cryptic promoter and non-coding cis-regulatory elements upon 5' deletions in NRXN1. We also observed recurrent intragenic deletions of ABCB11, encoding a transporter implicated in anti-psychotic response, in five treatment-resistant SCZ cases and showed that ABCB11 is specifically enriched in neurons forming mesocortical and mesolimbic dopaminergic projections. Our results indicate potential roles of sCNVs in SCZ risk.

G-quadruplexes associated with R-loops promote CTCF binding.

CTCF is a critical regulator of genome architecture and gene expression that binds thousands of sites on chromatin. CTCF genomic localization is controlled by the recognition of a DNA sequence motif and regulated by DNA modifications. However, CTCF does not bind to all its potential sites in all cell types, raising the question of whether the underlying chromatin structure can regulate CTCF occupancy. Here, we report that R-loops facilitate CTCF binding through the formation of associated G-quadruplex (G4) structures. R-loops and G4s co-localize with CTCF at many genomic regions in mouse embryonic stem cells and promote CTCF binding to its cognate DNA motif in vitro. R-loop attenuation reduces CTCF binding in vivo. Deletion of a specific G4-forming motif in a gene reduces CTCF binding and alters gene expression. Conversely, chemical stabilization of G4s results in CTCF gains and accompanying alterations in chromatin organization, suggesting a pivotal role for G4 structures in reinforcing long-range genome interactions through CTCF.

High-throughput Oligopaint screen identifies druggable 3D genome regulators.

The human genome functions as a three-dimensional chromatin polymer, driven by a complex collection of chromosome interactions1-3. Although the molecular rules governing these interactions are being quickly elucidated, relatively few proteins regulating this process have been identified. Here, to address this gap, we developed high-throughput DNA or RNA labelling with optimized Oligopaints (HiDRO)-an automated imaging pipeline that enables the quantitative measurement of chromatin interactions in single cells across thousands of samples. By screening the human druggable genome, we identified more than 300 factors that influence genome folding during interphase. Among these, 43 genes were validated as either increasing or decreasing interactions between topologically associating domains. Our findings show that genetic or chemical inhibition of the ubiquitous kinase GSK3A leads to increased long-range chromatin looping interactions in a genome-wide and cohesin-dependent manner. These results demonstrate the importance of GSK3A signalling in nuclear architecture and the use of HiDRO for identifying mechanisms of spatial genome organization.

Spatial and temporal organization of the genome: Current state and future aims of the 4D nucleome project.

The four-dimensional nucleome (4DN) consortium studies the architecture of the genome and the nucleus in space and time. We summarize progress by the consortium and highlight the development of technologies for (1) mapping genome folding and identifying roles of nuclear components and bodies, proteins, and RNA, (2) characterizing nuclear organization with time or single-cell resolution, and (3) imaging of nuclear organization. With these tools, the consortium has provided over 2,000 public datasets. Integrative computational models based on these data are starting to reveal connections between genome structure and function. We then present a forward-looking perspective and outline current aims to (1) delineate dynamics of nuclear architecture at different timescales, from minutes to weeks as cells differentiate, in populations and in single cells, (2) characterize cis-determinants and trans-modulators of genome organization, (3) test functional consequences of changes in cis- and trans-regulators, and (4) develop predictive models of genome structure and function.

3D genome, on repeat: Higher-order folding principles of the heterochromatinized repetitive genome.

Nearly half of the human genome is comprised of diverse repetitive sequences ranging from satellite repeats to retrotransposable elements. Such sequences are susceptible to stepwise expansions, duplications, inversions, and recombination events which can compromise genome function. In this review, we discuss the higher-order folding mechanisms of compartmentalization and loop extrusion and how they shape, and are shaped by, heterochromatin. Using primarily mammalian model systems, we contrast mechanisms governing H3K9me3-mediated heterochromatinization of the repetitive genome and highlight emerging links between repetitive elements and chromatin folding.

Cohesin-mediated loop anchors confine the locations of human replication origins.

DNA replication occurs through an intricately regulated series of molecular events and is fundamental for genome stability1,2. At present, it is unknown how the locations of replication origins are determined in the human genome. Here we dissect the role of topologically associating domains (TADs)3-6, subTADs7 and loops8 in the positioning of replication initiation zones (IZs). We stratify TADs and subTADs by the presence of corner-dots indicative of loops and the orientation of CTCF motifs. We find that high-efficiency, early replicating IZs localize to boundaries between adjacent corner-dot TADs anchored by high-density arrays of divergently and convergently oriented CTCF motifs. By contrast, low-efficiency IZs localize to weaker dotless boundaries. Following ablation of cohesin-mediated loop extrusion during G1, high-efficiency IZs become diffuse and delocalized at boundaries with complex CTCF motif orientations. Moreover, G1 knockdown of the cohesin unloading factor WAPL results in gained long-range loops and narrowed localization of IZs at the same boundaries. Finally, targeted deletion or insertion of specific boundaries causes local replication timing shifts consistent with IZ loss or gain, respectively. Our data support a model in which cohesin-mediated loop extrusion and stalling at a subset of genetically encoded TAD and subTAD boundaries is an essential determinant of the locations of replication origins in human S phase.

Cohesin-dependence of neuronal gene expression relates to chromatin loop length.

Cohesin and CTCF are major drivers of 3D genome organization, but their role in neurons is still emerging. Here, we show a prominent role for cohesin in the expression of genes that facilitate neuronal maturation and homeostasis. Unexpectedly, we observed two major classes of activity-regulated genes with distinct reliance on cohesin in mouse primary cortical neurons. Immediate early genes (IEGs) remained fully inducible by KCl and BDNF, and short-range enhancer-promoter contacts at the IEGs Fos formed robustly in the absence of cohesin. In contrast, cohesin was required for full expression of a subset of secondary response genes characterized by long-range chromatin contacts. Cohesin-dependence of constitutive neuronal genes with key functions in synaptic transmission and neurotransmitter signaling also scaled with chromatin loop length. Our data demonstrate that key genes required for the maturation and activation of primary cortical neurons depend on cohesin for their full expression, and that the degree to which these genes rely on cohesin scales with the genomic distance traversed by their chromatin contacts.

Chromatin domain alterations linked to 3D genome organization in a large cohort of schizophrenia and bipolar disorder brains.

Chromosomal organization, scaling from the 147-base pair (bp) nucleosome to megabase-ranging domains encompassing multiple transcriptional units, including heritability loci for psychiatric traits, remains largely unexplored in the human brain. In this study, we constructed promoter- and enhancer-enriched nucleosomal histone modification landscapes for adult prefrontal cortex from H3-lysine 27 acetylation and H3-lysine 4 trimethylation profiles, generated from 388 controls and 351 individuals diagnosed with schizophrenia (SCZ) or bipolar disorder (BD) (n = 739). We mapped thousands of cis-regulatory domains (CRDs), revealing fine-grained, 104-106-bp chromosomal organization, firmly integrated into Hi-C topologically associating domain stratification by open/repressive chromosomal environments and nuclear topography. Large clusters of hyper-acetylated CRDs were enriched for SCZ heritability, with prominent representation of regulatory sequences governing fetal development and glutamatergic neuron signaling. Therefore, SCZ and BD brains show coordinated dysregulation of risk-associated regulatory sequences assembled into kilobase- to megabase-scaling chromosomal domains.

Neuronal Yin Yang1 in the prefrontal cortex regulates transcriptional and behavioral responses to chronic stress in mice.

Although the synaptic alterations associated with the stress-related mood disorder major depression has been well-documented, the underlying transcriptional mechanisms remain poorly understood. Here, we perform complementary bulk nuclei- and single-nucleus transcriptome profiling and map locus-specific chromatin interactions in mouse neocortex to identify the cell type-specific transcriptional changes associated with stress-induced behavioral maladaptation. We find that cortical excitatory neurons, layer 2/3 neurons in particular, are vulnerable to chronic stress and acquire signatures of gene transcription and chromatin structure associated with reduced neuronal activity and expression of Yin Yang 1 (YY1). Selective ablation of YY1 in cortical excitatory neurons enhances stress sensitivity in both male and female mice and alters the expression of stress-associated genes following an abbreviated stress exposure. These findings demonstrate how chronic stress impacts transcription in cortical excitatory neurons and identify YY1 as a regulator of stress-induced maladaptive behavior in mice.

See(quence) and ye shall find: higher-order genome folding in intact single cells.

Payne et al. (2020) combine in situ imaging and ex situ sequencing via spatially resolved unique molecular barcodes to query higher-order genome folding patterns in intact single nuclei from mouse embryos and human fibroblasts.

Alteration of genome folding via contact domain boundary insertion.

Animal chromosomes are partitioned into contact domains. Pathogenic domain disruptions can result from chromosomal rearrangements or perturbation of architectural factors. However, such broad-scale alterations are insufficient to define the minimal requirements for domain formation. Moreover, to what extent domains can be engineered is just beginning to be explored. In an attempt to create contact domains, we inserted a 2-kb DNA sequence underlying a tissue-invariant domain boundary-containing a CTCF-binding site (CBS) and a transcription start site (TSS)-into 16 ectopic loci across 11 chromosomes, and characterized its architectural impact. Depending on local constraints, this fragment variably formed new domains, partitioned existing ones, altered compartmentalization and initiated contacts reflecting chromatin loop extrusion. Deletions of the CBS or the TSS individually or in combination within inserts revealed its distinct contributions to genome folding. Altogether, short DNA insertions can suffice to shape the spatial genome in a manner influenced by chromatin context.

3DeFDR: statistical methods for identifying cell type-specific looping interactions in 5C and Hi-C data.

An important unanswered question in chromatin biology is the extent to which long-range looping interactions change across developmental models, genetic perturbations, drug treatments, and disease states. Computational tools for rigorous assessment of cell type-specific loops across multiple biological conditions are needed. We present 3DeFDR, a simple and effective statistical tool for classifying dynamic loops across biological conditions from Chromosome-Conformation-Capture-Carbon-Copy (5C) and Hi-C data. Our work provides a statistical framework and open-source coding libraries for sensitive detection of cell type-specific loops in high-resolution 5C and Hi-C data from multiple cellular conditions.

Publisher Correction: Three-dimensional genome restructuring across timescales of activity-induced neuronal gene expression.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Three-dimensional genome restructuring across timescales of activity-induced neuronal gene expression.

Neuronal activation induces rapid transcription of immediate early genes (IEGs) and longer-term chromatin remodeling around secondary response genes (SRGs). Here, we use high-resolution chromosome-conformation-capture carbon-copy sequencing (5C-seq) to elucidate the extent to which long-range chromatin loops are altered during short- and long-term changes in neural activity. We find that more than 10% of loops surrounding select IEGs, SRGs, and synaptic genes are induced de novo during cortical neuron activation. IEGs Fos and Arc connect to activity-dependent enhancers via singular short-range loops that form within 20 min after stimulation, prior to peak messenger RNA levels. By contrast, the SRG Bdnf engages in both pre-existing and activity-inducible loops that form within 1-6 h. We also show that common single-nucleotide variants that are associated with autism and schizophrenia are colocalized with distinct classes of activity-dependent, looped enhancers. Our data link architectural complexity to transcriptional kinetics and reveal the rapid timescale by which higher-order chromatin architecture reconfigures during neuronal stimulation.